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We study the low frequency admittance of a small metallic island coupled to a gate electrode and 
to a massive reservoir via a multi channel tunnel junction. The ac current is caused by a slowly 
oscillating gate voltage. We focus on the regime of inelastic cotunneling in which the dissipation of 
energy (the real part of the admittance) is determined by two-electron tunneling with creation of 
electron-hole pairs on the island. We demonstrate that at finite temperatures but low frequencies 
the energy dissipation is ohmic whereas at zero temperature it is super-ohmic. We find that (i) the 
charge relaxation resistance (extracted from the real part of the admittance) is strongly temperature 
dependent, (ii) the imaginary and real parts of the admittance do not satisfy the Korringa-Shiba 
relation. At zero temperature the charge relaxation resistance vanishes in agreement with the recent 
zero temperature analysis [M. Filippone and C. Mora, Phys. Rev. B 86, 125311 (2012) and P. Dutt, 

T. L. Schmidt, C. Mora, and K. Le Hur, Phys. Rev. B 87, 155134 (2013)]. 
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I. INTRODUCTION 

During the last decades Coulomb blockade has become 
a powerful tool for observation of interaction and quan¬ 
tum effects in single electron devices m- This phe¬ 
nomenon is widely observed in low temperature electron 
transport through a single electron transistor. Another 
system which low temperature properties are affected by 
Coulomb blockade is a single electron box (SEB). It is 
schematically shown in Fig. Small metallic island is 
coupled capacitively to the gate electrode with the volt¬ 
age Ug. The number of electrons on the island is not 
conserved due to the tunneling in and out of an equilib¬ 
rium electron reservoir. A time dependent gate voltage 
Ug{t) generates ac current through the device. 

The equivalent electric circuit of a SEB (see Fig. 
is characterized by two capacitances. The gate capaci¬ 
tance Cg controls the external (induced) charge q on the 
island, q = CgUg. The total capacitance C determines 
the so-called charging energy Ec = e^/2C. It is the latter 
that is responsible for the Coulomb blockade effects. The 
tunnel junction is characterized by the dimensionless (in 
units e^/h) conductance g. Throughout the paper we use 
a standard assumption that the Thouless energy of the 
island is the largest energy scale in the problem. This 
allows us to work in a zero dimensional approximation 
neglecting spatial dependence of all quantities. 

Since there is no dc transport through the SEB, an es¬ 
sential dynamic characteristic becomes the admittance 
which characterizes the response of ac current I^ to 
the infinitely small ac part of the time dependent 
gate voltage Ug{t) = I/q + UujCosujt: Q{uj) = luj/U^. 
Long ago it was demonstrated that the admittance of 
SEB is affected by Coulomb blockade at low tempera¬ 
tures T ^ Ec [7]. However, since then the majority 
of works have addressed the so-called quantum capaci¬ 


tance: Ceff = dQjdUg, where Q is the average charge 
on the island, which determines the imaginary part of 
the admittance [8lll4]. Classically, at high temperatures 
T ^ Ec the effective capacitance coincides with Cg. As 
temperature decreases Ceff starts to deviate from Cg due 
to interaction and coherence effects. In seminal paper 
m it was suggested that the real and imaginary parts 
of the admittance in a SEB can be related in an universal 
way. After paper m a SEB admittance have attracted 
significant theoretical interest pRH2Q] . The admittance 
in the quasi-static regime was measured in a single chan¬ 
nel SEB constructed in 2D electron gas m- At present, 
there exists a number of measurements of admittances for 
different realizations of a SEB performed with the help 
of radio-frequency refiectometry [22H25] . 

The classical electrodynamics of a SEB suggests the 
following expression for the admittance at low frequen¬ 
cies, uj gEc'. 

g(u>) = -iOjCg + OJ^CgCR, (1) 

where R = h/{e^g) stands for the classical resistance of 
the tunnel junction. In Ref. the following general¬ 
ization of the classical result has been proposed for 



FIG. 1. (Color online) The set-up: a SEB subjected to a 
time-dependent gate voltage Ug{t) (left) and the equivalent 
electric circuit (right). 
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the quantum coherent SEB with C = Cg\ 

G{uj) = —icjCeff (2) 

where Rq was termed as charge relaxation resistance. 
Treating the Coulomb interaction within the Hartree- 
Fock approximation, the authors of Ref. m demon¬ 
strate that for single channel tunnel junction the charge 
relaxation resistance in Eq. § becomes universal, Rq = 
h/{2e‘^). The full quantum mechanical treatment of the 
charging energy in the case of single channel tunnel junc¬ 
tion demonstrates that at zero temperature Rq = h/ (2e^) 
{Rq = h/e^) for frequencies oj 6 (cu ^ S) [26]. Here 
S denotes the mean level spacing of single particle states 
inside the island of a SEB. Both results follow from two 
observations: (i) the effective low energy Hamiltonian 
of single channel SEB is of Fermi liquid type; (ii) the 
Korringa-Shiba relation for the response function 

iQ{ijj)luj holds within Fermi liquid low-energy description 
[26] . Recently, the analysis of Ref. [26| has been gener¬ 
alized to the case of a SEB with a weak ^ 1) multi 
channel tunnel junction and a large island, (5^0. It was 
found [29130] that at zero temperature the charge relax¬ 
ation resistance is inversely proportional to the number 
of channels in a tunnel junction and is independent of 
the external charge; Rq vanishes in the limit of infinite 
number of channels for any value of q. 

At finite temperatures the SEB with a multi channel 
tunnel junction in the limit of negligible mean level spac¬ 
ing, (5^0, has been analyzed in Ref. m- In particular, 
it was demonstrated that in the limit of weak tunneling, 
^ ^ 1, and near the charge degeneracy points the SEB 
admittance at low frequencies {uj <C ^max{|A|,T}) can 
be set down in the following form 

g{uj) = -icoCeS + LO^^CgTZq. (3) 

Here A denotes the electrostatic energy due to one excess 
electron on the SEB island. It depends on the external 
charge q and satisfy inequality |A| ^ Ec near a charge 
degeneracy point. The quantity Cg = dQjdUo stands for 
the renormalized gate capacitance which measures the 
response of the effective charge Q, introduced in Refs. 
[3M^ by one of us, to the static part of the gate volt¬ 
age. Contrary to the average charge Q on the island, 
the effective charge is expected to be integer quantized 
at zero temperature [32ll3^ . This implies that Cg van¬ 
ishes at T = 0 contrary to Ceff. The charge relaxation 
resistance in Eq. (§ is determined by the renormalized 
tunneling conductance g{T), IZq = h/{e^g{T)) ^ h/e^. 
The very same conductance g{T) determines the dc con¬ 
ductance of the single electron transistor (SET) under 
small bias between source and drain. We also note that 
Eq. © was proposed for the SEB with arbitrary relation 
between C and Cg. For the case of weak tunneling the 
treatment of Ref. m was restricted to the sequential 
tunneling approximation dressed by the renormalization 
due to virtual processes. The processes of inelastic cotun¬ 
neling [35| which dominate the dc transport through the 


SET at low temperatures T ^ Tin ^ were 

not taken into account. Therefore, the extrapolation of 
Eq. © down to the zero temperature and comparison 
with the result of Ref. [29130] were not possible. 

The real part of the admittance of a SEB with a 
multi channel tunnel junction in the regime of inelas¬ 
tic cotunneling has been studied in Ref. [7]. It was 
found that the real part of admittance is proportional 
to max{T^, This results implies the zero 

charge relaxation resistance at T = 0 in agreement with 
the result of Ref. [29130] extrapolated to limit of the infi¬ 
nite number of channels in the tunnel junction. However, 
the analysis of Ref. [7| has been restricted to Coulomb 
valleys, i.e. to integer values of q. 

In this paper we address the following question: how 
the zero temperature result for the charge relaxation re¬ 
sistance obtained in Refs. crosses over to the fi¬ 

nite temperature result of Ref. [31]? To answer this 
question we performed a detailed study of the admit¬ 
tance of the multi channel SEB near the charge neutrality 
points in the low temperature regime where the inelas¬ 
tic cotunneling processes dominate the dynamics. We 
found that the real part of admittance is proportional to 
max{T^, cj^l/A^. The charge relaxation resistance 
(extracted from Eq. Q) is strongly temperature depen¬ 
dent and small, Rq ^ (h/e^)(T/A)^ ^ hje^. In agree¬ 
ment with Refs. BSEn], we obtained that Rq is indepen¬ 
dent of g and vanishes at zero temperature. Our explicit 
results demonstrate strong violation of Korringa-Shiba 
relation for the response function iQ{(jj)luj and, conse¬ 
quently, support the non-Fermi liquid behavior of the 
multi channel SEB near the charge degeneracy points. 

The structure of the paper is as follows. In Sec. |TI| 
we introduce the Hamiltonian and Kubo formula for ad¬ 
mittance of a single electron box. The pseudofermion 
representation for the low energy Hamiltonian valid in 
the cotunneling regime is presented in Sec. m The re¬ 
sults of calculation of the admittance at low frequencies 
to the second order in the tunneling conductance g are 
given in Sec. El Finally, discussion of our results and 
conclusions are presented in Sec. [Vj The details of cal¬ 
culations are summarized in Appendix [^ We use units 
with h = e = 1 through out the paper except for the final 
results. 


II. FORMALISM 

A. Hamiltonian 

We start with the standard Hamiltonian describing 
Coulomb blockade in a SEB [36ll38] : 

H = Hi^HdEHcE Ht. 


(4) 
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Here Hi {H^) denotes free electron Hamiltonian in the 
lead (the island), 


where n^(co’) stands for the Fourier transform of the re¬ 
tarded polarization operator of electrons on the island: 


^ (5) 

a 

where operators a\ create an electron in the reser¬ 
voir (the island). Energies are counted from the 

chemical potential. The Hamiltonian 

Hc = Ec{nd-qf, = (6) 


n^(t) = i0(t)([nrf(t),nrf(O)]). (11) 

Here 0(t) denotes the Heaviside step function. In the 
quasi-static regime cj ^ 0, the polarization operator 
n^(u;) can be expanded in regular series in uj\ 

n^(cL;) = TTo + iuJiTi + O(cj^), (12) 

where both ttq and tti are real functions. We stress that 
we assume |ci;| ^6 throughout the paper. The static 
part TTo is fully determined by the average charge on the 
island Q\ 


takes into account the electrostatic energy due to the 
finite size of the island. 

The Hamiltonian 

-fft = E ^kaolda + h.C. (7) 

k,a 


describes tunneling of electrons between the island and 
the reservoir. In order to characterize the tunnel junc¬ 
tion, following Ref. m, we introduce Hermitian matrix: 


gaa’ = (27r)^ 




nl/2 


E4fe'^(4“^)w (8) 


TTo = 7^-C'eff - C, (13) 

where we remind Ceff = dQjdUQ with Q = {fid). This 
result holds by virtue of the Ward identity which re¬ 
lates the static polarization operator and the compress¬ 
ibility [39] . We note that on the classical level Ceff = Cg 
and TTo = 0. The classical electrodynamics result 0 im¬ 
plies TTi = 2'kC‘^ jg on the classical level. In what follows 
we discuss how quantum effects due to inelastic cotun- 
neling change naive classical expectations for ttq and tti. 


III. WEAK TUNNELING REGIME 


acting in the Hilbert space of the island’s states. Here 
the delta-functions are assumed to be smoothed on some 
intermediate scale between 6 and min{T, |cj|}. The ma¬ 
trix g allows one to define the number of open channels 
Nch and the effective channel conductance ^ch- 


(tr5)' 


trg2 


9ch — 


tr4 

tig ■ 


( 9 ) 


The dimensionless conductance g which characterizes the 
tunnel junction in classical electrodynamics is given as 

9 — ^ch-^ch¬ 
in the present paper we assume that Nch ^ 1 and 
^/^ch ^ 9ch ^ Although within these assumptions 
the classical conductance g can be still large, in what 
follows we restrict our consideration to the case 1 ^ ^ ^ 
1/^ch- We are interested in temperatures much smaller 
than the charging energy but much larger than the mean 
level spacing, Ec > T > (5. 


B. Admittance and polarization operator 


A. Projected Hamiltonian 

In what follows, contrary to Ref. [ 7 ], we confine our 
consideration to the vicinity of the one of charge de¬ 
generacy points, i.e. points where the external charge 
q = k ^ 1/2 where k is an integer. At these points the 
gap A = 2Ec{k + 1/2 — g) between the ground and the 
first excited state of the charging Hamiltonian He van¬ 
ishes. In the vicinity of the degeneracy point, the gap A is 
small in comparison with the charging energy, |A| ^ Ec. 
The processes of inelastic cotunneling becomes impor¬ 
tant at low temperatures T ^ |A|. At |A| ^ Ec, one 
can truncate the Hilbert space of electrons on the isolated 
island to two charging states characterized hy Q = k and 
Q = k ^ 1 [8]. The projected Hamiltonian acquires a 
form of 2 X 2 matrix acting in the isospin 1/2 space of 
these two charging states [8]: 


H = Hi+Hd + Ht + AS, + ^ + 


T 

4 ’ 


(14) 


where Hi^d are given by Eq. ([^, and 


The admittance of a SEB being the linear response of 
the ac current to the ac part of the time-dependent gate 
voltage, Ug{t) = Uo-\-Uuj coscjt, can be expressed as [31] 

g{uj) = -iujCg (1 + n^(u;)/c), (10) 


Et = ^ tkcxa\do,S^ - 1 - h.C.. (15) 

k,a 

Here E iS^ stand for standard spin 1/2 

operators. The Hamiltonian 0 is the Hamiltonian for 
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the A^ch channel Kondo model. The gap A between the 
charging states plays the role of an effective magnetic 
field. 

In the presence of ac component of the gate voltage, 
the energy detuning from the degeneracy point becomes 
time dependent: A(t) = A — {Cg/C)Uuj cosut. Then, 
as follows from Eq. ( p!4| ), the linear response of a SEB 
to ac gate voltage U^j is determined by the longitudinal 
dynamical isospin susceptibility: 

nf{t) = ie{t){[s,{t),SM])- (i6) 

In particular, the admittance is given as [ST]: 


= -ico^n% 


(17) 


We note that the average charge on the island can be 
written as Q = /c + 1/2 — {Sz)^ i.e. it is related to the 


isospin magnetization. As the consequence of Eq. (13), 


the dynamical isospin susceptibility should satisfy the re¬ 
lation nf (0) = CCefi/Cg. 


B. Pseudofermion effective action 


To deal with spin operators we employ the method 
of Abrikosov’s pseudofermion operators. Eollowing Ref. 
[40] . we introduce fermion operators t/Sq,, such that 

S = ^'0aO'a/3'0/3. (18) 

Here a = cr^, cr^} denotes standard Pauli matrices. 
To exclude redundant unphysical states (the states with 
> 1) wc add to the Hamiltonian H an arti¬ 
ficial chemical potential 77 for the pseudofermions. It is 
necessary to take the limit 77 ^ —00 at the end of any 
calculation. 

We remind that the physical partition function Z and 
correlation functions (O) can be found from the pseud¬ 
ofermion ones with the help of the following rules: 



In the case ^ 1 and 1 ^ ^ ^ after the 

integration over electrons in the reservoir and the island 
the Hamiltonian (14) can be transformed into the follow¬ 
ing imaginary-time effective action m- 


Q 

+ - / dTidT2a{Ti2)[lp{Ti)a-'lp{Ti)][lp{T2)a+ij{T2)]. 

( 20 ) 



FIG. 2. Feynman rules for the pseudofermion action (20). 


Here (7± = (cr^, ± and the kernel 


TT ' 


( 21 ) 


where ujn = 2TTTn is the bosonic Matsubara frequency. 
We note that the action similar to Eq. (20) has been first 


analyzed by Larkin and Melnikov in Ref. El- Effective 
action ( [^ corresponds to the XY Bose-Kondo model for 
the spin 1/2 [42143144] . 


The dynamical spin susceptibility (16) is determined 


by the pseudofermion dynamical spin susceptibility: 
n5,pf(r) = ^(7;[V5(r)cr^'^(T)][V;(0)cr^'^(0)]), 


( 22 ) 


where 7 / denotes the imaginary-time ordering and the 
average is taken with respect to the effective action ( 20 ). 


According to Eq. (19) we also need the expression for the 


physical partition function. It can be written as 




lim Z^ie 

?7^0 


-I3r] 






(23) 


where 0 cr(^) = stands for the exact 

imaginary-time pseudofermion Green’s function. 

As we discussed in the Introduction, in the Eermi liq¬ 
uid the imaginary and real parts of the dynamical spin 
susceptibility are related by the so-called Korringa-Shiba 
relation Taking into account that A in the effec¬ 

tive action (20) plays a role of magnetic field, Korring- 
Shiba relation for H^ should have the following universal 
form: 


Imnf(w) = 27rw[Renf(0)] , w ^ 0. (24) 


If this relation were correct, the low-frequency admit¬ 
tance Q{uj) = —iuj{Cg/C){CQ ^ iujC^Rq) would be char¬ 
acterized by the universal charge relaxation resistance, 
Rq = h/e^, similar to the single channel case. Here we 
introduce Co = (C/C^)Ceff. As we shall demonstrate 
below by direct calculation, Korringa-Shiba relation (24) 
does not hold for the effective action ( 20 ). 


IV. ADMITTANCE IN THE COTUNNELING 
REGIME 


The effective action (20) is suitable for the perturba¬ 


tion theory in ^ ^ 1 . Here we evaluate the imaginary 
part of the dynamical spin susceptibility n 5 (c(;) to the 
second order in g. The Eeynman rules for action (20) 
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(a) s + u + fl (b) 



FIG. 3. Feynman diagrams for the pseudofermion dynamical 
spin susceptibility in the hrst order in g. 


or hole remains on the island after each tunneling event. 
We remind that the SET conductance in the sequential 
tunneling approximation is also of the same Arrhenius- 
type form at low temperatures T ^ |A| ISEI]. As it is 
well-known, at low temperatures SET conductance has 
also a power-law (in T/|A|) contribution of the second 
order in g due to the processes of inelastic cotunneling 
[35]. As we shall demonstrate in the next section there 
is a similar contribution to the admittance. 


B. Perturbation theory in g: Inelastic cotunneling 


are shown in Eig.[^ The bare pseudofermion Matsubara 
Green’s function is given as follows: 


G — 


i£n-\-r] - crA/2‘ 


(25) 


Thus from Eqs. (22) and (23) in the zeroth order in g 
we find 




= 1 , 


= 2 cosh ——, 


nSf(*Wn)=0. (26) 


The processes of inelastic cotunneling becomes impor¬ 
tant at temperatures T Tin ^ |A|. However, these 
processes are of the second order in the tunneling con¬ 
ductance. Diagrams of the second order in g for the 
pseudofermion dynamical spin susceptibility are shown 
in Eig. [^ We remind that diagrams with pseudofermion 
loops vanish in the limit r] ^ — oo. Evaluation of the 
ten diagrams in Eig. [^yields the following result for the 
imaginary part of the low-frequency dynamical spin sus¬ 
ceptibility (see Appendix [A|) 


A. Perturbation theory in g: Sequential tunneling 


Before discussing the inelastic cotunneling regime (sec¬ 
ond order in g) we remind briefly the result of Ref. (31] 
for the pseudofermion dynamical spin susceptibility in 
the regime of sequential tunneling. In Eig. [^we present 
diagrams contributing to the pseudofermion dynamical 
spin susceptibility n 5 ^pf(ico’^) in the first order in g. Their 
evaluation demonstrates that ImHf (cj) suffers from sin¬ 
gularity at cj ^ 0: 




9 

Attuj sinh(;5A) 


(27) 


This unphysical divergence stems from non¬ 
commutativity of the limits cj ^ 0 and ^ ^ 0 in 
the structure of the ImH^ T<^)- Summing the ladder- 
type diagrams and taking into account finite (the lowest 
order in g) broadening of the pseudofermion Green’s 
function, we obtain the following expression m- 


nf(iad)(^) 


9 /3A 

47r sinh(;5A) 


gA ^ PA 

-luj + -— coth- 

27r 2 


-1 


(28) 

We note that the broadening of the singularity at cj = 0 
in Imn^pf(cc;) is determined by the sum of in- and out- 
rate of single electron tunneling [3] . At low temperatures 
T <C |A|, the result (28) implies ■(|w| < (7|A|/27r): 


* 5|A| 


(29) 


Imnf’(2)(w) = 


247rA4 


—V 

47r^ ) ’ 


1^1,r« |A|. 

(30) 


This expression dominates over the contribution (29) due 


to sequential tunneling at low temperatures, T ^ Tin- 
In the discussion above we do not consider renormal¬ 


ization of the effective action (20) between the ultra vio¬ 
let energy scale of the order of Ec and the infrared scale 
of the order of max{|A|,T} |4118145] . The renormalized 
effective action can be obtained from Eq. ( [^ by the fol¬ 
lowing substitutions p) g ^ g = Z‘^g 

and A A = Z^A, where the field renormalization fac¬ 
tor Z is given within one-loop approximation as im 


1 + 


In 


Er. 


1 - 1/2 


27r^ max{|A|,T}_ 


(31) 


We note that the pseudospin operator Sz renormalizes 
according to Sz Z‘^Sz [3T]. Then from Eq. (30) we 
find 


Iinnf>(2)(u;) = 


247rA^ 




(32) 


Remarkably, the field renormalization parameter Z drops 
out from Eq. (32) such that it coincides with Eq. (30) 


in spite of the renormalization. 

Taking into account one-loop renormalization of the 
effective action, one finds the following result for the av¬ 
erage charge on the island [8] : 


Such Arrhenius-type dependence, exp(—|A|/T), is char¬ 
acteristic of real processes in which an additional electron 




( 33 ) 
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factor. Then taking into account that (i) the relation 
dA/dA = Z‘^ holds within logarithmic accuracy and (ii) 
the renormalized conductance g corresponds to the en¬ 
ergy scale A, one can obtain the result (34). Alterna¬ 
tively, the appearance of the factor can be checked 
with the help of the expression for Q derived to the sec¬ 
ond order in g [33]. We emphasize that Eq. ( [M| ) implies 
that at T ^ Tin the effective capacitance becomes very 
different from Cg: Ceff = CgZ^gEj{27 t^\A\). 

Combining together Eqs. ( [M| ) and we obtain the 
following result for the admittance of a SEB (IcjI <^T 
Tin): 


r(\— ■ n ^ ^ ^ ^ ^ 

an - C, 


yoo) 


This is the main result of the present paper. Results 
( [ 3 ^ and (34) for the imaginary and real part of the 
dynamical isospin susceptibility implies strong violation 
of the Korringa-Shiba relation (24) at low temperatures 


when the processes of inelastic cotunneling dominate. 
Using the Korringa-Shiba relation, one overestimates er¬ 
roneously Imll^ (cj) by a large factor (A/T)^ ^ 1. The 
violation of the Korringa-Shiba relation signals that the 
Hamiltonian ( [l4| ) and, consequently, the effective action 
(20), involves non-Eermi liquid physics. 


V. DISCUSSION AND CONCLUSIONS 




FIG. 4. Second order in g diagrams for the pseudofermion 
dynamical spin susceptibility. 


In the case of low temperatures, T <C Tin, Eq. 
be simplified, and we obtain 


(33) can 


Renf’(i)(0) 


^ _ Z^g 

dA 47r^|A| 


(34) 


We note that the factor can be derived in the follow¬ 
ing way. Equation (31) determines the one-loop renor¬ 
malization group equation for the field renormalization 


According to Eq. ([^, the result ( [^ for the admit¬ 
tance implies the following result for the charge relax¬ 
ation resistance: 


Rc 


h TT^ / 


T 


3 V^^A. 


T<Ti, 


|A| 


ln(l/5)' 


(36) 


We emphasize that due to Coulomb interaction the 
charge relaxation resistance is strongly temperature and 
gate voltage dependent. Moreover, Rg depends on 
the tunneling conductance g although through the field 
renormalization factor only. Therefore, the charge relax¬ 
ation resistance depends in nontrivial way on the param¬ 
eters of a SEB in contrast to the zero temperature predic¬ 
tions of Refs. |29l3Qj and original ideas of Ref. [15]. Also 
the charge relaxation resistance is much smaller than the 
resistance quantum, Rq h/e^. At T = 0 the charge 
relaxation resistance vanishes, Rq = 0. This behavior is 
in agreement with the zero temperature result of Refs. 
|29l3Qj . Rq{T = 0) = (h/e^)/Ach, which implies zero 
charge relaxation resistance at T = 0 and Nch 00 . 

The real part of the admittance determines the energy 
dissipation rate of a SEB, W = CgEcReQ{uj)U^, which 
appears due to the time dependent periodic gate voltage 
modulations. The result (35) leads to the ohmic dissipa¬ 
tion at low frequencies: 


37r2A4 ^3^“’ 


|w| <C T < Tw 


(37) 
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FIG. 5. The processes of inelastic cotunneling. The interme¬ 
diate state has an additional electron (a) and hole (b). 


This result for the dissipation rate has the following 
physical explanation. Let us estimate the rate Fin for the 
two-electron process in which one electron with energy Si 
tunnels from the island into the reservoir and occupies 
the state with energy £4 whereas the other electron with 
energy £2 tunnels from the reservoir into the island occu¬ 
pying the state with energy 53 (see Fig. [^. For £i 7 ^ 53 
this process is inelastic and results in the electron-hole 
pair on the island at the end. It can go through two 
different intermediate states: with an additional electron 
and an additional hole on the island. The former costs 
the energy of the order of 2 Ec whereas the latter costs 
the energy of the order of |A|. In the considered case 
|A| ^ Ec, we can neglect the contribution due to the 
intermediate state with an additional electron. Provided 
such transition is accompanied by a periodic perturba¬ 
tion which supplies the energy cu to the final state, we 
can estimate the corresponding rate as follows 


II J dsjj fF{ei)[l - fries)] fF{£2) 

X [1 - /F(e4)](5(e3 +£4 - £1 - £2 - w). (38) 

Here we use the fact that typical electron or hole en¬ 
ergies in the integral in Eq. (38) are of the order of 
max{T, |ci;|} ^ |A|. Then for |c(;| ^ T ^ |A| we find that 
the frequency dependent part of the rate is estimated as 
f[^^ ^ g‘^T‘^uj/. In fact, this rate is similar to the rate 
derived for a SET biased by voltage in which case the role 
of uj is played by dc voltage [35]. The quadratic depen¬ 
dence of F-| f^ o n T is responsible for the factor in the 
expression (37) for the energy dissipation rate. Indeed, 
averaging T^^ over time-dependent gate voltage, taking 
contribution proportional to U‘^ and multiplying the re¬ 
sult by cj, one gets the estimate for the rate of energy 
dissipation, W ^ /dA^){CgU^^/CY, which, up 

to numerical factors, coincides with the result (37). At 


zero temperature the rate of inelastic cotunneling is given 


as Fin ^ g‘^ 00 ^!A?. This suggests the following estimate 
for the energy dissipation rate W ^ g^E^^CgU^/A^^ 
i.e. non-ohmic dissipation of energy at zero temperature. 
Using the result ( [3^ with T = 0, we obtain the following 
expression for the energy dissipation rate at zero temper¬ 
ature: 


W = 


Z^l^^2rj2 

127r4A4 9 


r«; |w| lAI. 


(39) 


It is worthwhile to discuss the result of Ref. [7] which 
is complementary to the result of the present paper. 
In Ref. [7] the real part of admittance was analyzed 
in the regime of inelastic cotunneling but precisely at 
q = i.e. at Coulomb valleys. It was found that 

Re0(ci;) (X max{T^,We emphasize that 
the result of Ref. [7] is quite unexpected. For example, 
the contribution due to inelastic cotunneling is propor¬ 
tional to T^/A^ near the charge neutrality points and to 
at Coulomb valleys [35]. Therefore, on the basis 
of our result one could expect that at Coulomb valleys 
the energy dissipation rate is given by Eqs. (37) and 


(39) with A substituted by Ec- However, this naive ar¬ 


gument leads to overestimation of the energy dissipation 
rate by large factor E‘1/ max{T^, cj^}. The result of Ref. 
[7] comes from particular cancellation of terms propor¬ 
tional to jE^ and bo^jE^ in Re0(co’) (see comment 
at the end of Appendix 0- „ 

We mention that our result ( j^ is at odds with the 
expression © proposed by us in Ref. m- Since the ef¬ 
fective charge Q is expected to be robustly integer quan¬ 
tized at zero temperature [32143^ . the renormalized gate 
capacitance Cg is exponentially small at T ^ |A| [33] . 
The renormalized conductance is known to be propor¬ 
tional to the temperature squared, g{T) ^ in 

the regime of the inelastic cotunneling [35133] . Therefore, 
Eq.(35) suggests the exponential suppression of the en¬ 


ergy dissipation at T ^ |A| contrary to the result (37). 
Thus for g ^ 1 the expression © works within the se¬ 
quential tunneling approximation dressed by renormal¬ 
ization due to virtual processes only. 

The effective action (20) predicts zero value of g un¬ 


der the renormalization in the infra red. The Nch chan¬ 
nel Kondo model (14) has the unstable fixed point at 
finite value of ^ ^ 1/A^ch [8146] . Therefore, our 

results obtained within the effective action ( 20 ) are ap¬ 


plicable for the Hamiltonian ( |14[ ) while g ^ g^. As fol¬ 
lows from Eq. ( jMj ), this condition implies that our re¬ 
sults hold not too close to the charge degeneracy point, 
|A| > igEc/'K'^) eye£>i-F‘^/gb}- Since g^ ~ 1/A^ch, the 
scale igEc/n^) exp(— tt^/^*) becomes extremely small al- 
ready for not too very large values of Nch- Comparing 
our result ( [3^ with the zero temperature result of Refs. 
[ 2 ^ , we find that for the case of finite number of chan¬ 
nels the charge relaxation resistance is given by Eq. (36) 
for temperatures T ^ \A\/^/N^. 

To summarize, we have studied the low frequency ad¬ 
mittance of a multi channel SEB under a slowly oscil¬ 
lating gate voltage. Eocusing on the regime of inelastic 
























cotunneling, we have calculated the admittance G(uj) (see 
Eq. ([^) near the charge degeneracy points. We found 
the following: 

(i) At finite temperatures but low frequencies, Tin ^ 
T ^ \uj\, the energy dissipation rate (determined by the 
real part of the admittance) is ohmic and scales as the 
temperature squared, see Eq. (37), in agreement with 
qualitative arguments. 

(ii) At zero temperature the energy dissipation rate 
is super-ohmic, ^ see Eq. (39), in agreement with 
qualitative estimates. 

(iii) The imaginary and real parts of the response func¬ 
tion iQ{uj)luj do not satisfy Korringa-Shiba relation. This 
supports the non-Eermi liquid behavior of the model near 
the charge degeneracy points. 

(iv) The charge relaxation resistance Rq is 
strongly temperature dependent and small, 
Rq ^ (h/e^)(T/A)^ ^ h/e^. It vanishes at T = 0 
in agreement with the recent zero temperature analysis 
of Refs. |29l3Qj . 

(v) The relation between the real part of admittance 
and the effective charge Q conjected by us in Ref. m 
does not hold beyond the sequential tunneling approx¬ 
imation dressed by renormalization due to virtual pro¬ 
cesses. 

Einally, we mention that our result (35) for the ad¬ 
mittance can be tested in a single electron box with 
small metallic island via radio-frequency reflectometry 
measurements mES]. Also we mention that following 
approach of Refs. our results can be extended to 

non-equlibrium conditions, e.g. different temperatures of 
the island and the reservoir l49l. 
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Appendix A: Computation of the polarization 
operator: Diagrams of the second order in g 

In this appendix we present details of computation of 
the polarization operator within the second order pertur¬ 
bation theory in g. There are contributions from the ten 
diagrams shown in Eig. The task is simplified consid¬ 
erably by the fact that we only need the imaginary part 
of the retarded polarization operator. Each diagram con¬ 
sists of six pseudofermion Green’s function lines and two 
interaction lines. Thus each diagram involves the sum¬ 
mation over three internal energies: fermionic 6 and two 
bosonic ones ft, ft'. As usual the fermionic sum is eas¬ 
ily undertaken with the help of the following identity: 

f{e) = (47ri)“^ f dstdinh{s/2T)f{s) where the con¬ 
tour of integration circles around all the poles of tanh. 


1. The diagram I 


The contribution from the diagram I to the polariza¬ 
tion operator can be written as 

cr=± 

xG(j(^is T —(jiis T iff)G—a-iis T ift -]- iuJn) 

xG(j{i£ T ift T ift ^G(j{i£ T ift T ift -t- icj^), (Al) 


where we introduce the kernel a{ift) = ga{ift)/A = 
^|r^|/(47r). Evaluating the sum over the fermionic energy 
5 and performing analytic continuation, iuj^ cj + iO, 
we obtain 


Imll 


s,pf 


(w) 


T 

Auj 




+ E 


/g - f-a 

2(jJ 


Im 




1 2 


(A2) 


Here = JpAv + Acr/2), /; = 

and /f(^) = 1/[1 + exp(5/T)] denotes the Eermi-Dirac 
distribution function. The functions are re¬ 

tarded function corresponding to the following Matsub- 
ara function 


n 


[a{ift + iuj) — a{ift)]^ 
{iuj)^{ift -f Aa)'^ 


(A3) 


2. The diagrams Ila and Ilb 


The contribution from the diagram Ila to the polariza¬ 
tion operator can be written as 

cr = ± 

xG(j{i£ T i^n)G‘^^{i£ T iff)G(j{i£ T ift T ift ). (-'^d) 

The contribution from the diagram lib can be found 
from the expression above by reverting the sign of uJn- 
{—iujn)• Evaluating the sum over 
the fermionic energy £, combining two contributions to¬ 
gether, and performing analytic continuation, we find 

= -Oa E 

cr=± 

+ f E T + f'-a) (A5) 


Here we introduce the following function 
Y = tW 


(A6) 
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Strictly speaking the sum defining is divergent. The 
summation is truncated at Q = flmax ^ Ec which is 
the model cut-off. The function T^r is therefore cut-off 
dependent. It is also important to note that when eval¬ 
uating summation over boson frequencies Vt' and coming 
across divergent expressions we assume symmetric limits 
~^max < < ^max- Tmucatcd Symmetric sums allows 

us to repeatedly shift a summation variable safely. 


3. The diagrams Ilia and Illb 


The contribution from the diagram Ilia to the polar¬ 
ization operator can be written as 

cr=± £,0,0' 

y^G(j{i£ T iujji)G—(j{i£ T i^^G—(jiis T 


TABLE I. The first few integrals Ik- 
h = +47r^T2)/3 

/2 =a)^(w^ + 47r2T2)/6 
h = -a)(3w^ + 207r^T2a)^ + 327r'‘T^)/30 
h = + 247r'‘T^)/15 


continuation, we find 


cr=± 


cjcr 


f-a 


cr = ± 


U) 




1 2 




cr = ± 


Auj 


(AlO) 


5. The diagram V 


The contribution from the diagram Illb can be found 
from the expression above by reverting the sign of 

Evaluating the sum over 
the fermionic energy 5, combining two contributions to¬ 
gether, and performing analytic continuation, we find 


The contribution from the diagram V to the polariza¬ 
tion operator can be written as 

ni^pf^(*W„) = -L y] cr^ y] a{in)a{in')Gl{ie) 

cr=± s,Q,Q' 

xG‘^(^i£ -|- iujji)G—(j{i£ -f- i^^^G—(jiis ff (All) 


im n 


fa — f-o 


s,pf 


cr=± 


2 uj 


Im 




1 2 


+ y] (A8) 


4. The diagrams IVa, IVb, IVc, and IVd 


The contribution from the diagram IVa to the polar¬ 
ization operator can be written as 


Evaluating the sum over the fermionic energy 5 and per¬ 
forming analytic continuation, we find 


Im n 


=- E 

cr=± 






a=zt 


2uJ(T 


f' 4- f' 

-T yy Im 


cr = ± 


Scj 


(A12) 


=-L y^(-cr 2 ) Y a{i9,)a{i9.’)Gl{ie) 

cr=± e,Q,0' 

xG(j{f£ H“ 'i^^n)G—(jiis A iCt^^G A 

xG—(j{f£ A A iujji)' 

The contribution from the other three diagrams can 
be found from the expression above by reverting 
the sign of uJn and the sign of a and the summa- 
tion sign: 11; 

T-r(2),IVc/. X ’ ■n-(2),IVa/. ’ x , 

ns,pf (^^n,cr) = n^pf [iuJn,-cr), and 

^s,pf (^^n,cr) = n;pf {-^uJn,(T). Evaluating 

the sum over the fermionic energy 5, combining all 
four contributions together, and performing analytic 


6. The result for Imnf^pf(ct;) in the 2d order in g 


Combining all contributions, Eqs. (A2), (A5), (A8), 
(AlO), and (A12), together, we find 


Im = -J2{K + f-.) 

a 

-2d A Y 


cr = ± 


UJCF 


_rj. Y ^ Im A:^’(24) 


UJ 


+2y AAz^im 


cr = ± 




(A13) 
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The following comment is in order here. Alternatively, 
one can compute the admittance by means of the current- 
current correlation function [ 7 ] . The latter consists of two 
operators which are of the first and seco nd or der in g [50] . 
Then the first and second lines in Eq. (A14) comes from 


the renormalization of the operator of the first order in 
g whereas the third and forth lines correspond to the 
contributions from the operator of the second order in g. 
Using Eq. (26), we find 


Iinnf>(2)(w) = -2tanh^ E 


cr=± 


cosh(/3A/2) 


5a E 

(A14) 


cr=± 


- W Im 

^ 12 w ^ . 


This expression is convenient for analysis at |u;|, T ^ | A|. 
Conv erting the sum over Matsubara frequencies in Eq. 
(A3) into the integral and performing analytic continua¬ 
tion, icj ^ cj + iO, we obtain 


Im7^5(1,i)(^)= J 


X ImQ!^(5 -h uj). 


(A15) 


Here = coth(5/2T), d^(cj) = —d^{uj) = —iguj/^An)^ 
and D^{e) = [5 + Act + i0]“^. Since lmD^{s) = 
— 7 r 6 {e + Act), the imaginary part of is expo¬ 

nentially small, ^ exp(—|A|/T). We note that 
is defined by the divergent Matsubara sum. It should be 


understood as a finite sum truncated at U = Un 


Er. 


Then the analytical continuation is possible. Moreover 
the imaginary part of is Ec independent. 

The terms proportional to are responsible for 

the renormalization of the first order perturb ative result 
(27). Therefore, only the last line in Eq. dAlil con¬ 
tributes to Imnfin the regime |cj|,T <C |A|. 


Again con verting the sun over Matsubara frequencies 
in Eq. (A3) into the integral and performing analytic 
continuation, iuj ^ uj iO, we find 


ImAT. 


CX) 

^’(2.1) = J - Be) Re -D^(e) Im a^{e) 


X Ima«(£ + a;) = (f^ (^16) 

^^TT/ ^7ru;2(Acr)^ ^ ^ 


Here we perform expansion in series in 1/A. The func¬ 
tions Ik are defined as follows 


oo 

h= J dee>^{e + uj){Be+u:-Be). (A17) 

— OO 

We list several first functions Ik (for k = 1, 2, 3,4) in Ta¬ 
ble jl} We note that the very same functions Ik determine 
the imaginary part of ’ T<^)- 


^ / g fc(-l)"-Vfc 

\47r/ 7rci;^(Acr)^+^ 


(A18) 


Using the functions Ik from the Table from Eq. (A14) 
we obtain the result (30) in the main text. 

Einally, we note that for the case of q = k, i.e. at 
the Coulomb peaks, one cannot adopt the pseudofermion 
technique and has to take into account all charging states. 
Then, the contribution to Re0(cc;) due to the in elastic 
cotunneling is given by the last line of Eq. (A14) multi¬ 
plied by a factor of 2 and with the following substitutions: 
A ^ En and 


(A19) 

uj uj \ 2Ec{2Ec — ooa) J 

It is due to this additional term the contribution of the 
order of 1 /E^ cancels and the admittance becomes pro¬ 
portional to 1/E’c H- 


^ G. Schon and A. Zaikin, “Quantum coherent effects, phase 
transitions, and the dissipative dynamics of ultra small 
tunnel junctions,” Phys. Rep. 198, 237 (1990) 

^ Z. Phys. B: Condens. Matter 85, 317 (1991, special issue on 
single charge tunneling, ed. by H. Grabert and H. Horner.). 
^ H. Grabert and M. H. Devoret, eds.. Single Charge Tun¬ 
neling (Plenum, New York, 1992). 

^ Y. Blanter and M. Biittiker, “Shot noise in mesoscopic con¬ 
ductors,” Phys. Rep. 336, 1 (2000) 

^ 1. Aleiner, P. Brouwer, and L. Glazman, “Quantum effects 
in Goulomb blockade,” Phys. Rep. 358, 309 (2002) 

® L. 1. Glazman and M. Pustilnik, “Low-temperature trans¬ 
port through a quantum dot,” in New Direetions in 
Mesoseopie Physies (Towards to Nanoseienee)^ edited by 


R. Fazio, V. F. Gantmakher, and Y. Imry (Kluwer, Dor¬ 
drecht, 2003). 

^ Yu. V. Nazarov, “Dissipation in a Goulomb capacitor,” 
Sov. J. Low Temp. Phys. 16, 422 (1990). 

® K. A. Matveev, “Quantum fluctuations of the charge of 
a metal particle under the Goulomb blockade conditions,” 
Sov. Phys. JETP 72, 892 (1991). 

^ H. Grabert, “Rounding of the Goulomb staircase by the 
tunneling conductance,” Physica B 194, 1011 (1994). 

H. Grabert, “Gharge fluctuations in the single-electron 
box: Perturbation expansion in the tunneling conduc¬ 
tance,” Phys. Rev. B 50, 17364 (1994) 

K. A. Matveev, “Goulomb blockade at almost perfect 
transmission,” Phys. Rev. B 51, 1743 (1995) 






















11 


X. Wang and H. Grabert, “Coulomb charging at large con¬ 
duction,” Phys. Rev. B 53, 12621 (1996) 

G. Goppert, H. Grabert, N. V. Prokof’ev, and B. V. Svis- 
tunov, “Effect of tunneling conductance on the Coulomb 
staircase,” Phys. Rev. Lett. 81, 2324 (1998) 

I. S. Beloborodov, A. V. Andreev, and A. I. Larkin, “Two- 
loop approximation in the Coulomb blockade problem,” 
Phys. Rev. B 68, 024204 (2003) 

M. Biittiker, H. Thomas, and A. Pretre, “Mesoscopic ca¬ 
pacitors,” Phys. Lett. A 180, 364 (1993) 

M. Biittiker and A. M. Martin, “Charge relaxation and 
dephasing in Coulomb-coupled conductors,” Phys. Rev. B 
61,2737 (2000) 

S. E. Nigg, R. Lopez, and M. Biittiker, “Mesoscopic charge 
relaxation,” Phys. Rev. Lett. 97, 206804 (2006) 

S. E. Nigg and M. Biittiker, “Quantum to classical tran¬ 
sition of the charge relaxation resistance of a mesoscopic 
capacitor,” Phys. Rev. B 77, 085312 (2008) 

Z. Ringel, Y. Imry, and O. Entin-Wohlman, “Delayed 
currents and interaction effects in mesoscopic capacitors,” 
Phys. Rev. B 78, 165304 (2008) 

H. C. Park and K.-H. Ahn, “Admittance and noise in 
an electrically driven nanostructure: Interplay between 
quantum coherence and statistics,” Phys. Rev. Lett. 101, 
116804 (2008) 

J. Gabelli, G. Eeve, J.-M. Berroir, B. Plagais, A. Cavanna, 
B. Etienne, Y. Jin, and D. C. Glattli, “Violation of Kirch- 
hoff’s laws for a coherent RC circuit,” Science 313, 499 
(2006) 

E. Persson, C. M. Wilson, M. Sandberg, G. Johansson, 
and P. Delsing, “Excess dissipation in a dingle-electron 
box: The Sisyphus resistance,” Nano Lett. 10, 953-957 
( 2010 ) 

G. Giccarelli and A. J. Eerguson, “Impedance of the single¬ 
electron transistor at radio-frequencies,” New J. Phys. 13, 
093015 (2011) 

S. J. Ghorley, J. Wabnig, Z. V. Penfold-Eitch, K. D. Peters- 
son, J. Erake, G. G. Smith, and M. R. Buitelaar, “Measur¬ 
ing the complex admittance of a carbon nanotube double 
quantum dot,” Phys. Rev. Lett. 108, 036802 (2012) 

J. G. Erake, S. Kano, G. Giccarelli, J. Griffiths, 
M. Sakamoto, T. Teranishi, Y. Majima, G. G. Smith, 
and M. R. Buitelaar, “Radio-frequency capacitance spec¬ 
troscopy of metallic nanoparticles,” Sci. Rep. 5, 10858 
(2015) 

G. Mora and K. Le Hur, “Universal resistances of the 
quantum resistance-capacitance circuit,” Nat. Phys. 6, 697 
( 2010 ). 

J. Korringa, “Nuclear magnetic relaxation and resonance 
line shift in metals,” Physica 16, 601 (1950). 

H. Shiba, “The Korringa relation for the impurity nu¬ 
clear spin-lattice relaxation in dilute Kondo alloys,” Prog. 
Theor. Phys. 54, 967 (1975). 

M. Eilippone and G. Mora, “Eermi liquid approach to the 
quantum RC circuit: Renormalization group analysis of 
the Anderson and Goulomb blockade models,” Phys. Rev. 
B 86, 125311 (2012) 

P. Dutt, T. L. Schmidt, G. Mora, and K. Le Hur, “Strongly 
correlated dynamics in multichannel quantum RC cir¬ 


cuits,” Phys. Rev. B 87, 155134 (2013) 

Ya. 1. Rodionov, 1. S. Burmistrov, and A. S. loselevich, 
“Gharge relaxation resistance in the Goulomb blockade 
problem,” Phys. Rev. B 80, 035332 (2009) 

1. S. Burmistrov and A. M. M. Pruisken, “Goulomb block¬ 
ade and superuniversality of the 0 angle,” Phys. Rev. Lett. 
101, 056801 (2008) 

1. S. Burmistrov and A. M. M. Pruisken, “Macroscopic 
charge quantization in single-electron devices,” Phys. Rev. 
B 81, 085428 (2010) 

A. G. Semenov, “On the macroscopic quantization in 
mesoscopic rings and single-electron devices,” (2013), 
arXiv: 1307.6615 

D. V. Averin and Yu. V. Nazarov, “Virtual electron dif¬ 
fusion during quantum tunneling of the electric charge,” 
Phys. Rev. Lett. 65, 2446 (1990) 

E. Mezei, “Theory of electron tunneling via real interme¬ 
diate states,” Phys. Rev. B 4, 3775 (1971) 

R. 1. Shekhter, “Zero anomalies in the resistance of a tunnel 
junction containing metallic inclusions in the oxide layer,” 
Sov. Phys. JETP 36, 747 (1973). 

1. O. Kulik and R. 1. Shekhter, “Kinetic phenomena and 
charge discreteness effects in granulated media,” Sov. Phys. 
JETP 41, 308 (1975). 

A. A. Abrikosov, L. P. Gorkov, and 1. E. Dzyaloshinski, 
Methods of Quantum Field Theory in Statistieal Physies 
(Prentice-Hall, Englewood Gliffs, NJ, 1963). 

A. A. Abrikosov, “Electron scattering on magnetic impu¬ 
rities in metals and anomalous resistivity effects,” Physics. 
Physique. Eizika. 2, 5 (1965). 

A. 1. Larkin and V. 1. Melnikov, “Magnetic impurities in an 
almost magnetic metal,” Sov. Phys. JETP 34, 656 (1972). 

S. Sachdev and J. Ye, “Gapless spin-fluid ground state in 
a random quantum Heisenberg magnet,” Phys. Rev. Lett. 
70, 3339 (1993) 

L. Zhu and Q. Si, “Gritical local-moment fluctuations in 
the Bose-Eermi Kondo model,” Phys. Rev. B 66, 024426 
( 2002 ) 

G. Zarand and E. Dernier, “Quantum phase transitions in 
the Bose-Eermi Kondo model,” Phys. Rev. B 66, 024427 
( 2002 ) 

H. Schoeller and G. Schon, “Mesoscopic quantum trans¬ 
port: Resonant tunneling in the presence of a strong 
Goulomb interaction,” Phys. Rev. B 50, 18436 (1994) 

G. Zarand, G. T. Zimanyi, and P. Wilhelm, “Two- 
channel versus infinite-channel Kondo models for the 
single-electron transistor,” Phys. Rev. B 62, 8137 (2000) 
Ya. 1. Rodionov, 1. S. Burmistrov, and N. M. 
Ghtchelkatchev, “Relaxation dynamics of the electron dis¬ 
tribution in the Goulomb-blockade problem,” Phys. Rev. 
B 82, 155317 (2010) 

Ya. 1. Rodionov and 1. S. Burmistrov, “Out-of-equilibrium 
admittance of single electron box under strong Goulomb 
blockade,” JETP Lett. 92, 696 (2010). 

Ya. 1. Rodionov and 1. S. Burmistrov, unpublished. 

Eshel Ben-Jacob, Emil Mottola, and Gerd Schon, “Quan¬ 
tum shot noise in tunnel junctions,” Phys. Rev. Lett. 51, 
2064 (1983). 


